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ABSTRACT 

We present a detailed analysis of a H 2 -bearing metal-rich sub-damped Lyman — a system at 
Zabs = 0.10115 towards the radio-loud quasar J0441—4313, at a projected separation of ~7.6 
kpc from a star-forming galaxy. The H 2 , Cl and Nal absorption are much stronger in the 
redder of the two components seen in the Hubble Space Telescope / Cosmic Origins Spectro¬ 
graph spectrum. The best single component fit to the strong H 2 component gives log ATfC) 
= 16.61 ± 0.05. However, possible hidden saturation in the medium resolution spectrum can 
allow for log ATH 2 ) to be as high as 18.9. The rotational excitation temperature of H 2 in 
this component is 133^22 K. Photoionization models suggest 30-80% of the total 7V(Hl) is 
associated with the strong H 2 component, that has a density 100 cm -3 and is subject to a 
radiation field that is ^ 0.5 times the Galactic mean field. The Very Large Baseline Array 1.4 
GHz continuum image of the radio source contains only 27% of the arcsecond scale emission. 
Using a previously published spectrum, no 21-cm absorption is found to be associated with 
the strong H 2 component. This suggests that either the A'(H 1 ) associated with this component 
is ^ 50% of the total iV(Hl) or the gas covering factor is SC 0.27. This is consistent with 
the results of the photoionization model that uses UV radiation due to stars in the associated 
galaxy. The 21-cm absorption previously reported from the weaker H 2 component suggests a 
spin temperature of SC 90 K, at odds with the weakness of H 2 , CI and Na I absorption in this 
component. From the inferred physical and chemical conditions, we suggest that the gas may 
be tracing a recent metal-rich outflow from the host-galaxy. 

Key words: galaxies: quasar: absorption line — galaxies: ISM — quasar: individual: 
J0441—4313 


1 INTRODUCTION 

Damped Lyman — a systems (DLAs) and sub-DLAs are by def¬ 
inition absorbers with neutral hydrogen column densit y, N( H I ) A 
2xl0 2 °cm -2 and ^ lxl0 19 cm~ 2 respectively (see lWolfe et al.l 
2005, for a review). These high Hi colum n density absorbers 
trace the bulk of the H I at 2 < 2 < 3 dPeroux et al.l 120051 : 
Noterdaeme et al. 2009bl. 120121) . and can in p rinciple contribute 
significantly to the global star formation rate |Wolfe et alj |2003l : 


* Based on observations made with the NASA/ESA Hubble Space Tele¬ 
scope, obtained from the data archive at the Space Telescope Science Insti¬ 
tute, which is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555; data obtained from 
the ESO Science Archive Facility under request number rduttal29485; data 
obtained with VLBA (Prgrm. ID: BD187, PI: Dutta) 
f E-mail: rdutta@iucaa.emet.in 


Srianand et al]2005l) . While they have been conjectured to be orig¬ 
inating from gas associated with high-* galaxies/proto-galaxies, 
their direct connection with galaxies is not yet well understood. 
The link between DLAs/sub-DLAs and galaxies can be established 
by directly detecting the host galaxies and/or showing that the pre¬ 
vailing physical conditions in the absorbing gas are consistent with 
those seen in a typical galactic interstellar medium (ISM). 


Due to difficulties associated with detecting high -2 galax¬ 
ies that are close to bright QSO sightlines, our understanding of 
the physical conditions in DLAs/sub-DLAs at 2 > 1.8 is pri¬ 
marily based on optical absorption line studies of low-ionization 
metal transitions a nd, in a few cases, H 2 , HD and CO molecu¬ 
lar transitions (e.g. Varshalovichet ai]|200ll: iLedoux et alJl2003l: 


iNoterdaeme et al.|[200 8a b. 2009a; Srianan d et alj 20081) . The ob 


served high-J excitations of H 2 at high -2 are consistent with 
a strong ultraviolet (UV) field suggesting in-situ star formation, 
while the typical temperatures and densities of the H 2 components 
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are 153 ± 78 K and 10—200 cm -3 respectively ISrianand et al.l 
1 20051) . The gas producing H 2 absorption pr obably traces diffuse 
molecular gas in the form of compact clouds (Balashe v et al.l201lh 
and contains only a small fraction of the Hi measured using the 
DLA profile dSrianand et al.ll2ol2l) . 

The HI 21-cm absorption line provides a complementary way 
to probe the physical conditions in the HI gas. If detected it can 
be used to investigate: (1) the thermal state of Hi gas, as 21-cm 
optical depth is a good tracer of the gas kinetic temp erature being 
invers ely proportional to the spin temperature dtCulkami & Heilesl 
mn>; (2) the parsec-scale structure in the absorbing gas via mil- 
liarcsecond (mas) scale spectroscopy (Sria nand et alj2013h : (3) the 
magnetic field in the cold neutral medium (CNM) using Zeeman 
splitting dHeiles & TrolandfeoQdl) : (4) the filling factor of cold gas 
in the ISM. Unfo rtunately, 21-cm detections from DLAs at high -2 
are very rare (see lSrianand et al1l2012l : lKanekar et al .1120 1 4[ for re¬ 
cent compilations). This lack of 21-cm detection is attributed to the 
gas being warm and/or to partial coverage of the background radio 
sources by the absorbing clouds. 

The effect of covering factor and the size of absorbing clouds 
can in principle be quantified by searching for 21-cm absorption 
in DLAs/sub-DLAs with H 2 detections. In our Galactic high lati¬ 
tude sightlines, a correspondence between H 2 and 21-cm absorp¬ 
tion is well established when log /V(H 2 ) ^ 15.8 (see IRov et al.l 
120061) . Based on this study, we know that the gas kinetic temper¬ 
ature measured from the rotational level populations of H 2 , T 01 , 
traces the Hi spin-temperature, T s . Therefore, in that case the lack 
of 21-cm absorption from a gas showing low rotational temperature 
of H 2 cannot be naturally attributed to the gas being warm. How¬ 
ever, high -2 DLAs with H 2 detection towards radio-bright QSOs 
(with flux density in excess of 100 mjy at the redshifted 21-cm fre¬ 
quencies), at appropriate 2 a bs (enabling 21-cm absorption search) 
are rare. Even the five cases where this is possible, do n ot follo w the 
relationship between 21-cm and H 2 absorption (see lSrianand et all 
120121 ). This lack of correspondence may mean that only a small 
fraction (i.e., ^ 10%) of the total HI is associated with the H 2 com¬ 
ponents or that the H 2 clouds are too small (i.e., < 15 pc) to cover 
even the unresolved radio components at mas scales. To resolve 
this, mas scale very long baseline interferometry (VLBI) imaging 
of the background radio source is e ssen tial to determine the cover¬ 
ing factor (see for example lGupta et aPl2012h . A good understand¬ 
ing of these issues is important to interpret the results of on-going 
and future 21-cm absorption line surveys, as well as surveys of H 2 
in DLAs/sub-DLAs. 

Low -2 DLAs and sub-DLAs are ideal targets to address the 
above issues. This is because: (1) for a given angular scale one will 
be probing a much smaller physical scale, where t he c overing fac¬ 
tor issues may have minimum effects (see lCurran et alii20051) : (2) 
in case of 21-cm detections one will be able to carry out VLBI 
spec troscopic observations to get spatial scales in the absorbing HI 
gas l lBorthakur et al.1120101 : ISrianand et aPl2013h : (3) direct asso¬ 
ciation of the absorbers with the galaxies responsible for the ab¬ 
sorption is possible. The last one will allow us to relate the in¬ 
ferred spatial scales in the absorbing gas to its location with re- 
spect to the galaxy (i.e., in the star forming regions or in the halo). 
ISrianand et alj ( 120130 have demonstrated the presence of structures 
in the 21-cm optical depth on parsec scales in the ISM of a low -2 
galaxy. In such cases of the background radio source being struc¬ 
tured, the spatial variation of the Hi gas opacity can be studied 
using VLBI spectroscopy. 

UV spectroscopic surveys with the Hubble Space Telescope 
(HST), are identifying significant numbers of low -2 DLAs and sub- 


DLAs jRao et alj 1201 it [Meiring et al J 1201 ll : iBattisti et al.l l2012h . 
Thanks to the Cosmic Origins Spectrograph (COS) on board the 
HST, there have been t hree recent cases of H 2 detection in low- 
2 (i.e., 2 < 0.6) DLAs JCrighton et alj|2013l : lOliveira et alj|2014l : 
Srianand et al]2014l) . Recently by a careful search in 27 low -2 (i.e., 
2 < 1) DLAs and sub-DLAs towards 26 QSOs, using HST/COS 
archival spectra. iMuzahid et al.ld2014l) have reported 7 new H 2 de¬ 
tections. With the largest sample to date (10) of H 2 detections at 
low- 2 , they find > 2 times higher H 2 incidence rate at low -2 com¬ 
pared to high- 2 . The H 2 components in their sample have typical 
temperatures of 133 ± 55 K. From photoionization models and the 
lack of high-J excitations, the authors have shown that the prevail¬ 
ing radiation field is much weaker than the Galactic UV radiation 
field, in contrast with the results at high- 2 . Moreover, the large im¬ 
pact parameters of the host-galaxies (> 15 kpc) for majority of their 
systems strongly suggest that the H 2 bearing gas is not related to 
star-forming disks but probably stems from extended halo gas. 

Among the low -2 sample of H 2 absorbers of IMuzahid et all 
HH, two are towards radio-loud QSOs, which in principle facil¬ 
itates detailed study of the cold gas in these absorbers. However, 
the redshift of one of these systems falls in the frequency range af¬ 
fected by radio frequency interference at available radio telescopes. 
The other system is the subject of detailed study in this paper. This 
paper is organized as follows. In Section [2] we give details of the 
observations and data reduction process. The physical properties 
of the system as derived from the absorption lines are discussed in 
Section[3] The properties of the associated galaxy and their compar¬ 
ison with that of the absorber are discussed in Section[4] Section[5] 
explains the results from our photoionization modelling of the ab¬ 
sorber. In Section [6] we focus on the Very Long Baseline Array 
(VLBA) mas scale imaging of the background radio source at 1.4 
GHz and the previously reported 21-cm absorption from the cold 
gas using the Australia Telescope Compact Array (ATCA) . Lastly, 
we summarize our results in Section [7] Throughout this paper we 
use a flat cosmology with Ho = 70 km s -1 Mpc -1 and = 0.27. 


2 OBSERVATIONS AND DATA REDUCTION 

The QSO J044117.3-431343 (2 em = 0.59378; also known as 
PKS0439—433 and henceforth J0441—4313) was observed us¬ 
ing COS during HST cycle-19, under program ID: 12536 (PI: 
V. Kulkarni). These observations consist of G130M far-UV grat¬ 
ing integrations (5.3 ks) at a medium resolution of R ~ 18000 
(full width at half maximum, FWHM ~ 18 kms -1 ). The data 
were retrieved from the HST archive and reduced using the CAL- 
COS pipeline software. The individual reduced xld files were first 
flux calibrated. Then the alignment and co-addition of the sepa¬ 
rate exposures were carried out using the software developed by 
Danforth et alj j20ldR The exposures were weighted by the inte¬ 
gration time while co-adding the flux calibrated data. The final co¬ 
added spectrum covers the observed wavelength range 1131 —1435 
A and has a signal-to-noise ratio, S/N ~ 16—20 per resolution el¬ 
ement. Since each COS resolution element is sampled by six raw 
pixels, we binned the spectrum by 3 pixels, which further improves 
the S/N per pixel. All measurements and analyses were, subse¬ 
quently, performed on the binned data. Measurements are, however, 
found to be fairly independent of the binning. Continuum normal- 


1 http://casa.colorado.edu/danforth/science/cos/costooIs.html 
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ization was done by fitting the line free regions with a smooth lower 
order polynomial. 

The COS wavelength calibr ation is known to be uncertain 
at the level of 10—15 km s -1 dSavage et all 1201 it [Meiring et al.l 
120131) . Regions of the spectrum that are recorded near the edges of 
the detector segment are likely_to have an erroneous wavelength so¬ 
lution (see for e.g. lMuzahid et al.l2014h . However, in the case of the 
present spectrum, the velocity offsets at different wavelengths as 
estimated from the difference in the observed and the expected line 
centroids of the numerous H 2 absorption lines are < ±5 kms -1 . 
Moreover, the Galactic absorption lines are found to be aligned 
within < ±5 kms -1 . The line spre ad function (LSF) of the COS 
spectrograph is not a Gaussian. iGhavamian et al.l J20091 ) provided 
a characterization of the non-Gaussian LSF for COS, which was 
subsequently updated by Kris 1 J201 ll ). We adopt this updated LSF 
for our Voigt profile fitting analysis. Interpolated LSF at the line 
centre was convolved with the model Voigt profile while fitting an 


In addition, we use the pipeline calibrated spectrum (R ~ 
45000, FWHM - 6.6 kms -1 ) of QSO J0441-4313 obtained us¬ 
ing the Ultraviolet Echelle Spectrograph (UVES) on the Very Large 
Telescope (VLT), available in the European Southern Observatory 
(ESO) archive 0 The spectra corresponding to different exposures 
covering our regions of interest, after applying barycentric correc¬ 
tion, were brought to their vacuum values using the formula given 
in lEdlenl Il966h . For the co-addition, we interpolated the individual 
spectra and their errors to a common wavelength array, and then 
computed the weighted mean using the weights estimated from the 
error in each pixel. The final spectrum covers the observed wave¬ 
length range 3732—6837 A with S/N ~ 40—50 per resolution ele¬ 
ment. 

Furthermore, we observed the background radio source with 
the VLB A using the 21-cm receiver band for 1 hr on 2014 August 
26. The total bandwidth was 256 MHz in dual polarization (eight 
32 MHz baseband channel pairs). Each baseband channel was split 
into 128 spectral points. Two-bit sampling and a correlator integra¬ 
tion time of 2 s were used. The observations were carried out in 
nodding-style phase referencing with a cycle of ~5 min, i.e., ~3.5 
min on source and ~ 1.5 min on the phase calibrator (J0440—4333). 
A strong fringe finder/bandpass calibrator (DA 193) was also ob¬ 
served at the beginning for 4—5 min. The target source was ob¬ 
served for ~35 min, split into scans at different hour angles to 
improve the UV-coverage. Data were calibrate d and imaged using 
Astronomical Image Processi ng System (AIPS; Greisen 200 ^) in a 
standard way (see for example fMomiian et al .1200 2; Srianan d et al.l 

|2oT3) . 


3 ABSORPTION LINE PROPERTIES OF THE SUB-DLA 

The sub-DLA at * a bs = 0.10115 towards J0441—4313 with 
log Af(Hl) =19.63 ± 0.08 and log iV(H 2 ) = 16.61 ± 0.05 
flMuzahid et al.ll2014j ) presents a very interesting case for under¬ 
standing the cold gas that is present around galaxies. This system 
was originally selected as a DLA candidate based on strong Mg II 
and F e II absorption in the Faint Object Spectrograp h (FOS) spec¬ 
trum dPetitiean et al J 1 9961) . We note that lChen et al.l d2005h derived 
log N (H I) = 19.85 ± 0.10 for this sub-DLA using the Space Tele¬ 
scope Imaging Spectrograph (STIS) G140L spectrum, consistent 

2 http://www.ast.cam.ac.uk/~rfc/vpfit.html 

3 http://archive.eso.org/wdb/wdb/adp/phase3_spectral/form 


within 2a of our measurement. In this section, we discuss the prop¬ 
erties of the sub-DLA as probed by the absorption lines (metals and 
H 2 ) detected in the COS and UVES spectra. 


3.1 Analysis of metal lines 


We detect the metal transitions C I, C II, N I, N II, O I, Si II, Si III, 
P II, S II, Arl. Fell and FeUI in the COS spectrum spread over 
a velocity range of ~150 kms -1 . Some of these absorption pro¬ 
files are shown in the left panel of Fig. Q] The ClI, N II, Ol, Sill, 
Si III and Fe III lines are saturated/heavily blended and hence not 
used in our analysis. We do not cover the C II* AA1335.7, 1335.6 
transitions, and the C II* A1037 transition is heavily blended. C I* 
absorption is absent and we place a 3 a upper limit on its column 
density using the strongest transition that is free from any blend, 
i.e., A1279.9. Further, Mg II AA1239, 1240 and Mnll A1197 ab¬ 
sorption are not detected and we use these to place 3a upper limits 
on N (Mg II) and JV(Mn II) respectively. From the Voigt profile fits 
to the unblended and most likely unsaturated metal lines detected 
in the COS spectrum, we infer two main absorbing clumps of gas at 
2 abs = 0.10094 & 0.10119 (henceforth Component 1 and 2 respec¬ 
tively), ~ 68 kms -1 apart (see left panel of Fig. |T}. We assume 
that all the neutral and singly ionized atoms are physically asso¬ 
ciated with the same gas cloud. Hence while fitting, the redshift 
and b parameter for each absorption component are tied to be the 
same for each of the ions, i.e., only turbulent broadening is consid¬ 
ered. Additionally, we rejected fits with more than two components 
since, although having similar , they have larger values of AICC 
(Akaike informatio n criteria 1 Akaike 1974) corrected for the finite 
sampl e size (AICC: ISugiuralll978h as given in equation 4 of lKind 
([ 2 OI l])) and higher errors in the parameters. 

The medium resolution of the COS spectrum is not sufficient 
to resolve all the different components in the absorbing gas, in par¬ 
ticular the narrow ones. In the present case, since Call and Nal 
absorption lines have been detected in the high resolution UVES 
spectrum dRichter et alJhoiH) . we can get a realistic idea of the 
component structure. We find that ten and eight components give 
the best fit to the Call and Nal lines respectively (see Table [7] 
and right panel of Fig.jT}. The column densities that we obtained 
are consistent within errors with those reported by iRichter et al .1 
l l201ll) . Table [2] lists the column density measurements or 3 a up¬ 
per limits in case of non-detections of the various ions. In the case 
of Ca II and Na I, to facilitate comparison between Components 1 
and 2 seen in the COS spectrum, we sum the column densities in 
components (a) to (d) and (e) to (j) respectively. 

Both the Arl lines (A1066 and A1048) are not detected in 
Component 1. Absorption from Component 2 at the expected wave¬ 
length ranges of Ar I is present. The Ar I A1066 line is blended with 
some other absorption, and the Arl A1048 line centroid is shifted 
by 1 pixel with respect to the expected peak absorption from other 
metal lines (see Fig. |TJ. It is not clear whether this shift is caused 
by wavelength scale uncertainties or not. Regardless, it is clear that 
the Ar I strength we find is lower than what we expect if there is no 
relative depletion between Ar and S and if ionization corrections 
are ignored. Later we will discuss this issue in detail. 

Since N( HI) cannot be decomposed into components, the av¬ 
erage abundance can be estimated as, [X/H] = log ( N (X N)//V(H I)) 
— log (A(X)/V(Hl ))0 + IC, where N(X N) is the sum of the 
column densities of all the components of the dominant ion N of 
element X, and IC is the average ionization correction. Through¬ 
out this paper we use the solar elemental abundances as given in 
lAsplund et al.1 d2009l) . In Table [2] we list the average metallicities 
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without any ionization corrections for those species which are ex¬ 
pected to be the dominant ions of the respective elements in the 
neutral phase. IC for the S abundance, determined from SII, is 
found to be negligible compared to the errors as discussed in Sec¬ 
tion [5] Here, IC = log fki/fsn, where fm and fsn are the ion¬ 
ization fractions of HI and S II respectively. Moreover, since S is 
known not to be affected by dust depletion, we can take the aver¬ 
age metallicity of the absorber as that of the S abundance, which 
is found to be super-solar ([S/H] = 0.28 ± 0.08). Additionally, the 
abundance of P, most likely a non-refractory element, is found to 
be super-solar. This is higher than the prev iously reported value 
of log Z = —0.20 ± 0.30 by IChen et aO (1 2005 ), based on the 
Fe II column density measured using the low-resolution FOS spec¬ 
trum, log A r (Hl) = 19.85 ± 0.10, and the mean measured value of 
IV(Zn Il)/A r (Fe II) in the DLA population. 

When we consider the total column densities, we find [Ar/S] = 
— 1.04 ± 0.14. The ratio [Ar/S] is ^ —1.66 and —0.85 in Compo¬ 
nents 1 and 2 respectively. This is much lower than the mean value 
of —0.4 found bv IZafar et alj ( j2014h for high-z DLAs. The large 
Ar depletion we find towards our high metallicity system with low 
N(H I) is in line with the mild correlati ons of the [Ar/a] ratio with 
metallicity and A(H I) as noted bv IZafar et al.l d2014l) . We wish to 
point out that such large Ar I depletions are seen in t he case of High 
Velocity Clouds (HVCs) in the Milky Way (see lRichter et al.l20011) . 
We discuss this issue in detail in Section [5] 

Fe II is the only iron co-production species clearly detected in 
our spectrum. In the absence of a Zn II column density measure¬ 
ment it is difficult to interpret the Fe abundance, since differences 
between Fe and any other a element may reflect either Fe deple¬ 
tion or nucleosynthetic origin. Here we proceed with the assump¬ 
tion that the lower abundance of Fe is a reflection of dust depletion. 
Hence, the dust depletion obtained is, [Fe/S] = —0.49 ± 0.04, the 
column density of dust in Fe is, log ./V(Fedust) = 15.24 ± 0.01, 
and the dust-to-gas ratio is, k = 1.29 ± 0.2 l O The inferred Fe de¬ 
pletion is much less than what is typically seen in the cold gas in 
the Galactic ISM. From Fig. 9 of Iweltv et al.l Jl996h . we notice 
that in the Galactic ISM, A r (Nal)/A r (CaIl) > 1 for the observed 
values of .(V(Nal) in the present system. On the contrary, we find 
that ATGall) is higher than ATNal) in the present case. This could 
mean Ca depletion is not as high as typically seen —3 dex) in 
the Galactic ISM. Such a picture is consistent with relatively small 
depletion we infer for Fe from the observed [Fe/S]. As noted pre¬ 
viously, we do not detect Mg II AA1239, 1240 and Mnll A1197. 
However, the 3a limits we obtain are not stringent enough to pro¬ 
vide further insights into the chemical history of this system. 

High-s DLAs with such high value s of m etallicities, 
^(Feaust) and k tend to show II 9 absorption dLedoux et al.ll2003l : 
IPetitiean et alj|2006l : iNoterdaeme et al.|l2008al) . H 2 molecular ab¬ 
sorption is indeed detected from this s ystem at 3 a b s = 0.10115 from 
J = 0, 1.2,3 levels bv lMuzahid et al.l i2014l) . We also notice a con¬ 
sistent albeit weaker H 2 absorption at 3 a bs = 0.10091. From equiv¬ 
alent width measurements of MgII AA2796, 2803, Fell A2600 
given i nlChurchilll (1200 lh using the FOS spectrum, and the results 
of lGuptai'et al.l i 2012lk we find that there is a high probability of de¬ 
tecting 21-cm absorption from this system. A tentative weak 21-cm 
absorption has been reported from this system at 2 a bs = 0.10097, 
however no 21 -cm absorption corresponding to_the stronger H 2 
component (i.e., Component 2) is detected 1 Kanekar et alJuOOlh . 


4 Here, lV(Fe duat ) = N(S)[1 - l 0 l Fe / s ]](Fe/S) 0 
and, k = 10l s / H l[l - 10[ Fe / s ]] 


Table 1 . Component-wise column densities (with errors shown in parenthe¬ 
ses) of Ca II and Na I detected from the sub-DLA at z a b a = 0.10115 towards 
J0441—4313 


Comp. 

z abs 

b 

(kms -1 ) 

log IV(Call) 
(cm -2 ) 

log 7V(Na I) 
(cm -2 ) 

a 

0.10084 (0.000001) 

1.55 (0.46) 

11.61 (0.04) 

11.26 (0.03) 

b 

0.10090 (0.000003) 

4.91 (1.03) 

11.63 (0.09) 

11.16(0.06) 

c 

0.10093 (0.000003) 

4.44(1.08) 

11.74(0.08) 

10.59 (0.25) 

d 

0.10100 (0.000005) 

8.80 (2.34) 

11.53 (0.08) 

10.46 (0.23) 

e 

0.10108 (0.000002) 

5.88 (0.93) 

11.48 (0.06) 

11.19(0.04) 

f 

0.10112(0.000002) 

1.55 (0.87) 

11.54 (0.06) 

— 

CT 

& 

0.10114 (0.000001) 

2.19(0.20) 

12.08 (0.03) 

12.04 (0.02) 

h 

0.10118 (0.000002) 

4.37(1.16) 

11.72 (0.05) 

11.33 (0.06) 

i 

0.10122 (0.000005) 

0.84 (0.80) 

11.37 (0.16) 

10.67 (0.16) 

j 

0.10125 (0.000007) 

5.35 (2.56) 

11.34 (0.13) 

— 


In Section [?7T2] we will discuss the implications of this 21-cm mea¬ 
surement in detail. 

We do not find any strong (i.e., > 3a) difference in the ra¬ 
tios N (Fe U)/N(S II), N(P II)/N(S II) and N (N I)/V(S II) between 
the two components. While most of the metal column densities in 
the two components match within ~0.2 dex, C I, Na I, Ar I and H 2 
absorption are much weaker in Component 1. Since these species 
(apart from Arl) can be ionized by photons with energy <11 eV, 
their column densities will be useful in constraining the ionizing 
radiation field, as explored in Section[5] 

From the column densities of N I, P II, S II, Call, Fell, we 
find that ~ 55—65% of the total metal column density is in Com¬ 
ponent 2. This implies that the fraction of /V(H I) in the strong H 2 
component is, /n(hi) ~ 0.55—0.65, if the metallicity is uniform 
across the components. On th e other hand, if A'(Nal) scales with 
/V(H I) as seen in our Galaxy dFerlet et alJll985l : IWakker &~M athis 
l 2000 l) . then we would expect ~80% of the total N (Hi) to be 
present in Component 2. Hence, from the observed distribution of 
the metal column densities across the two components, /n(hi) is 
expected to be ~ 0.55—0.8 in Component 2. If the known rel ations 
of A~ ( Nal) and ./V(Call) with A(Hl) in our Galaxy dFerlet et al.l 
1 19851 : IWakker & Mathis 200(J) were valid for the present system, 
we would expect log A r (Nal) = 11.49 and log A(Call) = 11.77. 
However, the observed column densities are ten times higher. In 
other words, ^10 times more Nal and Call per Hi are present in 
this sub-DLA than what is typically seen in the Galactic ISM. As 
both Nal and Call are not the dominant ions of the respective ele¬ 
ments in the HI phase, this may imply that the background ionizing 
field in this sub-DLA is weaker than the mean Galactic radiation 
field. We come back to this point in Section[5] 


3.2 Physical Conditions 

3.2.1 H 2 absorption 

Details of the H 2 absorption properties in the two components 
are given in Table [3] while the fits are shown in Fig. C2 of 
iMuzahid et al.l (120141) 7 We consider the AXH 2 ) obtained from VP- 
FIT for the weaker component as an upper limit as most of the tran¬ 
sitions are blended. Both the components are ~10 kms -1 shifted 
from the redshifts of the corresponding metal components derived 
using the COS spectrum. However, the stronger H 2 absorption 
arises within ~2 kms -1 of the strongest components of Call and 
Nal. The inferred excitation temperature is ^ 200 K in both the 
components, similar to what is expected from a CNM gas. How¬ 
ever, this gas is situated outside the visible optical disk of the can- 


© 0000 RAS, MNRAS 000, 000-000 




















































Cold parsec-scale gas 5 


Table 2. Component-wise column densities (with errors shown in parentheses) of metals detected from the sub-DLA at z a b s =0.10115 towards J0441—4313 


Ion (X N) 


log N(cm~ z ) 


[X n/H] c 

Component l a 

Component 2° 

Total 

Hi 

— 

— 

19.63 (0.08) 

— 

Cl 

13.31 (0.14) 

14.01 (0.04) 

14.09 (0.04) 

— 

Cl* 

+ 13.40 

+ 13.50 

+ 13.75 

— 

Nl 

14.64 (0.02) 

14.71 (0.02) 

14.98 (0.02) 

-0.48 (0.08) 

Nal 1 

11.60(0.04) 

12.18 (0.02) 

12.28 (0.02) 

— 

Mg 11 

+ 15.30 

+ 15.30 

+ 15.60 

+ 0.37 

PlI 

13.00(0.11) 

13.29 (0.07) 

13.47 (0.06) 

+0.43 (0.10) 

Sn 

14.65 (0.04) 

14.79 (0.03) 

15.03 (0.02) 

+0.28 (0.08) 

Arl 

+ 12.27 

13.22 (0.12) 

13.27 (0.14) 

-0.76 (0.16) 

Call 1 

12.24 (0.04) 

12.45 (0.03) 

12.66 (0.02) 

— 

Mn 11 

+ 12.80 

+ 12.90 

+ 13.2 

+ 0.14 

Fell 

14.49 (0.05) 

14.72 (0.03) 

14.92 (0.03) 

-0.21 (0.09) 


“ z aba = 0.10094 ± 0.000003 & b = 19.6 ± 1.5 kms" 1 
b 2abs = 0.10119 it 0.000003 & b = 23.8 ± 1.3 kins -1 
c Average metallicities without applying any ionization corrections 

1 The column densities in Components 1 & 2 are summed over components (a) to (d) & (e) to (j) respectively as listed in Table |T| 
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Figure 1. A selection of metal lines associated with the sub-DLA at z a bs = 0.10115 towards J0441—4313 detected in the COS (Left) and UVES (Right) 
spectrum. Best-fitting Voigt profiles are overplotted in red. The tick marks show the component positions. The errors in flux and residuals from the fit are 
shown at the bottom as green lines and blue histograms respectively. 


didate host-galaxy (see Fig. 1 of lChen et al.ll20()T If we consider 
the iV(H 2 ) of the weaker component as a measurement, the total 
Eh column density for this system is log IVfEh) = 16.64 ± 0.05, 
and hence the average molecular fraction is log f(Eh) = —2.69 
± 0.0® HD absorption is undetected and from the strongest un¬ 
blended transition of HD(J=0) (A 1042), we estimate a 3 u upper 

limit of log jV(HD,J=0) l< 13.90. _ 

As pointed out by iMuzahid et alj ( I20l4) . due to the moder¬ 
ate spectral resolution of the COS spectrum, N(H 2 ) measurements 
can be highly uncertain if the intrinsic 6 values are small. To test 
the uncertainties in the N{ H 2 ) measurement, we fitted the H 2 in 
Component 2 with multiple components with their z and b fixed 
to that of the Nal components arising within 30 kins -1 of z aba = 
0.10115. We find that the H 2 lines can be adequately fitted with at 


5 f(H 2 ) = 2JV(H 2 ) / (2JV(H 2 ) + 1V(HI)) 


most three of these components, and further components are super¬ 
fluous. Of these three components, the component corresponding to 
the strongest Na I with 6 = 2 km s _1 , contributes Js 99% to the total 
W(H 2 ). The resultant log IV(H 2 ) = 18.9 ±0.1 from this fit is much 
higher than that obtained from our single component fit with 6 as a 
free parameter. While the T 01 (127 ± 10 K) obtained from this fit is 
consistent with that obtained from our single component fit. How¬ 
ever, the multi-component fit gives higher AICC value (1675) as 
well as larger parameter errors compared to our single component 
fit. 

Note that our approximation of fixing the 6 of the H 2 lines as 
that of the Nal components is not realistic, especially if the ther¬ 
mal broadening is large. In addition, as the strongest Na I compo¬ 
nent contains almost all the N{ H 2 ), we can approximate the H 2 ab¬ 
sorber as a single cloud, and use the curve of growth (COG) to get a 
better handle on the b parameter of H 2 . For details_of the COG tech¬ 
nique to estimate the b of H 2 transitions see lSrianand et alj d2014t) . 
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We notice that the absorption lines of the J = 3 transitions from 
Component 2 which have a spread in the Xf values (A: rest wave¬ 
length, /: oscillator strength), show a wide range of rest-equivalent 
widths. Hence using a single cloud COG for H 2 ( J=3) transitions, 
we obtain b = 4—7 kms -1 for Component 2. This is of the or¬ 
der of one-third to one-half of the resolution of the COS spectrum. 
We then fit all the H 2 lines by fixing b to be 4—7 kms -1 in VP- 
FIT. Here we have assume d that the b values of d ifferent J levels 
are the same (see however [Noterdaeme et al.ll2007l) . These fits give 
log ATH 2 ) = 18.70—18.90 for Component 2, consistent with the 
multi-component fit with b fixed to that of the Nal components, 
and T 01 = 113 ± 10 K. However, the AICC for these fits are larger 
compared to the fit in which the & is a free parameter, indicating that 
the present data do not support a narrow b and a high iV(H 2 ). More¬ 
over, H 2 absorption from the higher J-levels (J ^ 4) as well as HD 
absorption, as expected for such large NQI 2 ), are not detected in 
the present spectrum. However, we note that a higher resolution 
spectrum is required to accurately estimate both these quantities. 
Here we proceed assuming log IVXH 2 ) = 16.61 ± 0.05 for H 2 in 
Component 2, though we explore the possibility of the gas having 
high N(H 2 ) using photoionization modelling in Section HU The 
important point to note here is that the T 01 does not seem to de¬ 
pend on the b parameter, and both the multi-component fit and the 
single cloud COG approach lead to T 01 consistent with our original 
estimate. 

We estimate a 3 a upper limit of log JV(H 2 , J=4) ^ 13.82, us¬ 
ing the strongest unblended line (A1044) available. This allows us 
to constrain the photoabsorption rate of Lyman- and Werner-band 
UV photons by H 2 in the J = 0 level (/3 0 ) from the equilibrium of 
the J = 4 level population. From equation (5) and values given in 
iNoterdaeme et al.|d2007L Section 4.2), we get fio ^ 2 x 10~ n s -1 , 
i.e., one-tenth of the Galactic rate or lower. From simple formation 
equilibrium of optically thin H 2 tlurall 19751) . we have the neutral 
hydrogen density, riH = 0 . 11 /?oW(H 2 )/(iiiV(Hl)), where R is the 
formation rate of H 2 . We use the ISM value of R (3 x 10 -17 s ” 1 
cm -3 ) scaled by the measured dust content in the system, to obtain 
the average density of the gas, riH ^ 50 cm~ 3 . Hence, from the 
H 2 formation, the average gas density is constrained to be less than 
~50 cm -3 , if the values of relevant parameters are similar to those 
measured in the Milky Way. 


3.2.2 Ionization and fine-structure excitation ofC° 

Assuming solar relative abundance of C and S, and that C II in the 
neutral phase can be traced by SlI, we can obtain WfCll). Un¬ 
der this assumption, we derive some of the physical conditions 
in Component 2 using standard techniques. From photoionization 
equilibrium between CI and CII and taking t he gas temperature as 
T 01 , we can estimate the electron density, n e dSrianand et al.l2005i . 
Equation (5)). Assuming the Galactic p hoto ionization rate for Cl 
(T = 2— 3.3 x 10~ 10 s~ 1 : |Peauignot & Aldrovandil d 1986h ) results 
in n e ~ 0.03—0.10 cm' 3 for Component 2. If C is depleted with 
respect to S, we expect the electron density to be higher than this 
value. On the other hand, if the radiation field strength is \H0 th of 
the Galactic mean field, as suggested from the lack of H 2 absorp¬ 
tion from J = 4 level, the inferred electron density will be lower. 
Taking njnn « 10~ 3 , typical of CNM gas dWolfire et alJI 1995h . 
leads to «,h ~ 30—100 cm' 3 in this gas for the above assumed T. 
Considering Component 1, the n e /r ratio is 3.6 times lower than 
that of Component 2. If we assume T to be similar for both the 
components, this will imply that the gas density in Component 1 is 
lower than that of Component 2. 


Table 3. Component-wise summary of H 2 absorption properties detected in 
th e su b-DLA at z a b s = 0.10115 towards J0441—4313 from lMuzahid et al] 

(m! 



Component 1 

Component 2 

■^abs 

0.10091 (0.000006) 

0.10115 (0.000001) 

log 1V(H2) (cm” 2 ) 

^ 15.51 (0.03) ° 

16.61 (0.05) 

b(kms _1 ) 

Toi(K) 

32.7 (2.6) 

178 t 37 6 

12.0 (0.5) 

1 72 +33 
—22 


“ Should be treated as an upper limit as this component is severely blended 
b Should be taken as indicative of the typical temperature expected in this 
component 



T (K) 


Figure 2. The ratio N(Cl*)/N(Cl) plotted as a function of temperature 
for different njj. The horizontal line shows the observed upper limit on 
1V(C l*)/N(C I), while the solid and dashed vertical lines are the inferred 
temperature and 1 -cr range respectively for the strong H 2 component. 


In addition, from CI and C I* equilibrium we estimate tih cor¬ 
respondin g to Co mponent 2, following the procedure and values 
given in lSrianand et aL ( 2005 1 Section 4.2). In Fig. [2] we plot the 
N(C I*)/iV(C I) ratio as a function of temperature for different nu- 
From the observed upper limit of this ratio and Tot estimated for 
Component 2, we get a limit on the density, nn ^38 cm -3 . While 
we have considered the Galactic UV pumping rate (Toi = 7.55 x 
10 _1 ° s- 1 ) here, we note that scaling the rate does not affect the 
result significantly, which is expected only when collisional excita¬ 
tions dominate. Therefore, the ionization and fine-structure excita¬ 
tion of CI are consistent with the gas density in Component 2 being 
few tens of H atoms cm -3 . Note that the N(C l*)/N(C I) ratio mea¬ 
sured for Component 1 is not stringent enough to place any useful 
constraints on the physical conditions in this gas. 
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4 CONNECTING GALAXY AND ABSORPTION 
PROPERTIES 


The sub-DLA towards J0441—4313 is associated with a spiral 
galaxy at an impact parameter of ~ 7.6 kpc from the QSO sightline 
l lPetitiean et al.ll 199a : IChen et al]|2005b . This is the lowest galaxy 
impact parameter at whic h Ha absorption has been detected at 
low-a (see Fig. 9 of IMuzahid et al J 12014f ). enabling the connec¬ 
tion between absorber and galaxy properties to be well-established. 
IChen et all d2005b have carried out spectroscopic analysis of the 
galaxy, and found it to be a star-forming galaxy with an oxygen 
abundance more than solar ([O/H] = 0.45 ± 0.15). They derived 
the sub-DLA metallicity as log Z = —0.20 ± 0.30 (see Section l+TV 
and hence reported a metallicity decrement along the galaxy disk. 
Using higher-resolution data, however, we have found the average 
metallicity in the sub-DLA to be twice solar. Therefore, there does 
not appear to be any considerable metallicity gradient (—0.02 ± 
0.17 dex kpc -1 ) along the disk of the galaxy. This is similar to the 
results of lPeroux et ID {mi, who do not find any significant de¬ 
crease between metallicity of galaxy from emission and metallicity 
of absorber from absorption, and even report the possibility of an 
increasing metallicity gradient in few cases. 


Using the luminosities of the galaxy's emission lines provided 
by IChen et all 1 20051) , and the abundance determination method 
i Izotov et ali i 20061) . we estimate [N/H] = —0.16 and hence 


given in 

[N/O] = —0.61 in the galaxy. Here we have assumed a typical elec¬ 
tron temperature of 10 4 K. The average [N/S] ratio for the sub-DLA 
is —0.76 ± 0.03 (ionization corrections are negligible compared 
to the errors). Note that this ratio is similar within 0.1 dex across 
the two components. The abundances of S and O, both being a- 
elements, track each other. Hence, the similarity of the [N/a] ratio 
in both the galaxy and the sub-DLA implies that the chemical en¬ 
richment history of both are similar. Thus, the absorber may be trac¬ 
ing the extended neutral gas disk of the galaxy. Indeed, rotational 
ve locity measurements of the galaxy and the sub-DLA performed 
by IChen et al.[ {2005i ) hints that that absorbing gas may be corotat¬ 
ing along with the optical disk at a galactocentric radius of ~ 13.6 
kpc (i.e., the deprojected separation of the sub-DLA along the stel¬ 
lar disk). Alternatively, the gas probed by the sub-DLA could be 
recently ejected from the galaxy and tracing a metal-rich galactic 
wind/outflow. 


Nitrogen can be of both primary and secondary origin, de¬ 
pending on whether the seed C and O are produced by the star 
itself during helium burning (primary), or whether they are from 
yields of earlier generations of stars and hence already present in 
the ISM from which the star formed (secondary). In nearby galax¬ 
ies, it has been observed that at lower metallicities, i.e., for [O/H] 
^ 0.4, the [N/O] ratio remains constant (primary N), while it rises 
steeply with increasing O abundance, i.e., for [O/H] > 0.7 (sec¬ 
ondary N). Therefore, the measured [N/a] in this sub-DLA along 
with the high metallicity places it in the secondary regime of N pro¬ 
duction. Hence, the absorbing gas must have undergone a substan¬ 
tial period of star formation. This reinforces our above hypothesis 
of either the galaxy having an extended neutral disc out to at least 
~13 kpc or the sub-DLA tracing the metal-enriched halo gas of the 
galaxy. We test both these scenarios through our photoionization 
models in Section[5] 

Another possible indication that the galaxy might have under¬ 
gone recent periods of outflow is the presence of a Lya system 
at z a bs = 0.10204, ~ +250 kms -1 from the sub-DLA (see Fig. 
A2 of lMuzahid et al.|[2014l) . with Lya, Si II, Si III and N II transi¬ 
tions spread over ~ 200 km s -1 . We measure log JV(Hl) = 14.93 ± 


0.44, log N(Si II) = 13.81+0.14, log N (Si III) = 13.99 ± 0.19, and 
log N (N II) = 14.81 ± 0.28 in this system. Based on the large metal 
column densities, the system could be metal-rich, with metallicity 
and [N/a] ratio close to solar. In that case, this cloud is likely to 
be tracing a recent metal-enriched outflow from the galaxy. How¬ 
ever, since only the Lya transition is covered for this system, the 
A’(Hl) measurement has large uncertainties and it is plausible that 
the system has higher N(H I) and lower metallicity. 


5 PHOTOIONIZATION MODELS 

We use the+hotoionization code CLOUDY (version 13.03; last de¬ 
scribed bv lFerland et alj2013l) to model the physical conditions and 
infer the chemical enrichment in this system. In all our models, the 
absorber is considered to be a plane-parallel slab of constant den¬ 
sity gas with the radiation field impinging on it from one side. Note 
that we do not aim to model all the absorbing ions detected from 
this system due to the multiple-component nature of the absorber 
and complexities associated with phase structure, depletion, and 
possible hidden saturation in the COS spectrum. We concentrate on 
H 2 , C I, C I*, and Na I absorption arising from Component 2, which 
are likely to trace the CNM and whose excitations will be governed 
by similar ioni zing radi ation. The ‘atom H2' command of CLOUDY, 
as described in lShaw et ali (I2005I) . is used in order to get an accu¬ 
rate H 2 equilibrium abundance. We consider three different inci¬ 
dent radiation fi elds: (i ) the me tagala ctic UV background of QSOs 
and galaxies lHaardt & Madaull2012l . herea fter HM12), (ii) the in¬ 
terste llar radiation field as in our Galaxy jBlack & van Dishoeckl 
11987b , scaled by a factor xuv, and (iii) a starburst galactic radia¬ 
tion field. Examples of the typical incident continuum spectra for 
different radiation fields used in our models are shown in Fig. [3] 
The cosmic microwave background radiation at z = 0.1 and a cos- 
mic ray ionization rate of log(rcfl) = —17.3 {williams et all 1998b 
are also included. For the observed N{ H I) , we find the ionization 
correction for the S abundance to be negligible (i.e., IC ^ — 0.02) 
compared to the errors, for a wide range of density and radiation 
fields considered. Hence, the average metallicity of the system can 
be represented by [S/H] = 0.28 ± 0.08, i.e Z ~ 2 Zq . 

In the following sections, we focus our efforts on modelling 
the strong H 2 component (Component 2) using CLOUDY simula¬ 
tions, that self-consistently compute the ion and molecular abun¬ 
dances along with the gas temperature. We run a grid of models 
by changing the density, ?ih (Sections l5. 1 l & [5~2l l. or the ionization 
parameter, U (Section l5.3t . and stopping the calculations when the 
7V(H2) in the model reaches the observed value of 10 16,61 cm -2 . 
Additionally, we discuss the effects of a high A r (H 2 ) value of 10 18,9 
cm -2 (see Section [3.2. Il l on the models in Section 15+1 We con¬ 
strain the density using the observed limit of N(C l*)/A r (C I), the 
N(Cl)/N(Nal) ratio, and by requiring the ./V(Hl) of the model to 
be less than the total measured ./V(Hl). Note that, as discussed in 
Section [XT1 fits to the metal lines in the COS spectrum with more 
than two components are not preferred. However, fitting the CI 
lines with multiple components having the z and b fixed as that 
of the Na I lines, allow for two times larger N(C I) than that ob¬ 
tained from our best fit. We allow for this uncertainty in the column 
density measurements while comparing with the model predictions. 
Since there are many parameters in the models, for simplicity we 
take the metallicity in this component to be the same as the average 
(2 Zq) and the H 2 formation rate, R, to be the same as that given by 
IJural {l975b for our Galaxy. Dust composition is assumed similar to 
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Figure 3. Examples of typical incident continuum spectra for different ra¬ 
diation fields used in our photoionization models. 

the Galaxy and dust depletion is taken as observed for this system. 
The results of the models are summarized in Table0 

5.1 Models with metagalactic UV background 

We first consider the scenario in which the absorber is a gas cloud 
irradiated with the HM12 background. From the total iV(Hl) con¬ 
straint, we find that the observed TVfFh) can be produced at ?ih ^ 
0.03 cm -3 , while from the upper limit of NIC 1*)/N(C I), we 
have riH < 50 cm -3 . For riH ~ 0.1 cm -3 , the N(C I)//V(NaI) 
ratio, N(Ci) and TV(Nal) predicted by the model is consistent 
with the observed values within the uncertainties. Moreover, the 
Tor predicted by this model is consistent with our estimation (see 
Table 0. Flence, just the extragalactic background can explain 
the observations provided the gas is at low densities. However, 
we note that this gas cloud is located just outside the optical 
disk of a star-forming galaxy (see Section 0. If we consider an 
incident spectrum similar to that of this galaxy and account for the 
measured dust extinction (see Section [53l for details), the expected 
radiation field near the absorber at 10 eV is ~10 times higher 
than the HM12 background. If the C I* absorption could be better 
constrained with higher resolution and S/N data or the C II* ab¬ 
sorption were covered, it would place additional constraint on the 
density, and hence allow us to test this model better. The fraction of 
7V(H I) associated with the H 2 gas in this case is, /n(hi) ~0.5. We 
estimate the size, L, of the H 2 absorbing cloud, assuming spherical 
geometry, as ~74pc, using L = / N (hi) N(Hi) / (f H i n H ). 


5.2 Models with the mean Galactic radiation field 

Next, we consider the case where the absorber is situated in the ex¬ 
tended neutral disc of a galaxy with radiation field similar to that 
of the Milky Way. The Galactic radiation field as given in ‘table 
ism’ of CLOUDY is the unextinguished local interstellar radiation 
field. However, most of the radiation field between 1 to 4 Rydbergs 
is found to be heavily absorbed by gas in the ISM. Hence, using 
the ‘extinguish’ command of CLOUDY we introduce photoelectric 
absorption by a slab of cold neutral gas (log iV(H I) ~ 20) to mimic 


typical Galactic ISM sightlines. This is a more appropriate way to 
model the CNM phase in the gas, since this will lead to a single 
neutral phase throughout the gas slab. Moreover, while the hydro¬ 
gen ionizing photons emitted from the galaxy are in principle re¬ 
moved, the species of our interest (H 2 , C I, Nal) are susceptible to 
ionizing photons with energy less than 1 Rydberg. Note that we add 
the HM12 radiation field to the Galactic radiation field. 

We find that when thermal equilibrium is assumed the T 01 es¬ 
timated by the models is too low (^30-40 K) compared to the ob¬ 
served T 01 (133l 3 2 K). The fact that thermal equilibrium models 
in CLOUDY produce very low T 01 , while using the radiat ion given 
by the ‘table ISM’ in CLOUDY, has been noted by [Srianand et ali 
d2014l) while modelling H 2 absorption in a low-3 DLA. They 
suggested that the low temperatures in the models could be due 
to CLOUDY not considering additional non-radiative heating pro¬ 
cesses. Additio nal heating in D LAs can c ome from cosmic-ray 
ionization (see lDutta et al.ll2ol4h . However. ISrianand et al.l ^20141) 
have shown that increasing the cosmic-ray ionization rate, while 
increasing the temperature, also causes the ion column densities 
to increase. Instead, they present H 2 formation at lower densities 
through enhanced formation rate on dust grains as a poss ible solu ¬ 
tion for increasing the gas temperature (see also lHabart et al.l201 it 
iLe Bourlot et al]2012h . 

Since the thermal equilibrium models are not able to 
self-consistently explain the observed temperature, we consider 
constant temperature models, i.e., we fix the temperature to be 
equal to the observed T 01 throughout the gas cloud. From these 
models using the iV(Hl) constraint and the 7V(C I*)//V(C I) limit, 
we find the H 2 can be produced for riH ~ 3—30 cm -3 and xuv ~ 
0.1—0.5. However, the model-predicted N(C I)/iV(NaI) ratio for 
the above solution does not match the observed ratio within the 
allowed uncertainties (see Table 0. Hence, while the Galactic 
radiation field attenuated by 10—50% can produce the observed 
N{ H 2 ), it is unlikely to produce the Cl and Nal absorption as 
measured for this gas, and a weaker field may be required. 


5.3 Models with starburst radiation field 

Next, we consider the absorber to be present in the halo around a 
galaxy. IChen et ahl J2005I) have measured properties of the galaxy 
believed to be associated with this sub-DLA (see Section0. Using 
their measurements, we construct a model with radiation field sim¬ 
ilar to that of the candidate host-galaxy. From the Ha extinction- 
corrected luminosity, we estimate a star formation rate of 0.53 
M@yr _1 dKennicutal 19981) . Subsequently, we use STARBURST9fQ 
to generate the spectrum of a continuously star-forming galaxy. The 
metagalactic HM12 UV background is added to the above radia¬ 
tion field. We proceed similarly as above to estimate the physical 
parameters of the system, except in this case we run the model on 
ionization parameter grids. As in Section [5~2l the T 01 estimated by 
this model when assuming thermal equilibrium is much lower than 
our measurement. Hence, we consider a model with the tempera¬ 
ture kept constant (i.e., equal to T 01 ). We find that such a model 
is consistent with our observations for —5.5 ^ log U ^ —5.0. This 
can be seen from Table0and Fig.0 where we have plotted N( H I), 
N(C I*)/A r (C I), 7V(Cl)/JV(NaI) and column densities of Cl and 
Na I obtained from the model as a function of ionization parameter. 
The N(C I) and JV(Na I) are over-predicted by the model by 0.3 dex 

6 http://www.stsci.edu/science/starburst99/docs/table-index.html 
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Figure 4. Results from photoionization model using the starburst radiation field for the strong H 2 absorbing component (Section l.73f The model-predicted 
IV(H I), JV(C l*)/N(C I), N(C l)/lV(Na I) and column densities of CI and Nal as a function of ionization parameter are shown in panels (a), (b), (c) and (d) 
respectively. The vertical shaded regions show the range of U in which the model is consistent with our observations. The horizontal shaded regions in (a), (b) 
and (c) show the allowed range in the values from observations. The symbols circle and square in (d) show the measured Y(C I) and ,Y( Na I), with the size of 
the symbols representing the allowed uncertainties in the measurements. 


and 0.5 dex respectively, which can be accounted for by depletion 
as seen in the Galactic CNM iWelly et all 19991 ). Over the range of 
U where the model is consistent with our observations allowing for 
uncertainties in measurements and depletion effects, /n(hi) is in 
the range ~ 0.3—0.8. 

We define the ionization parameter as, U = 
Qfle~ T /(AnCnHc), where Q is the rate at which ionizing 
photons are emitted by the galaxy, fl is the solid angle subtended 
by the galaxy at the absorber, r is the dust optical depth, and r 
is the distance between the galaxy and the illuminated face of 
the cloud. Q can be obtained by integrating the STARBURST99 
spectrum. Since we are estimating U using C I, Na I and H 2 
excitations, we estimate the dust optical depth at 10 eV (rio e rz)- 
From the extinction measurement of the galaxy, E(B — V) = 
0.22 dChen et all 1200511 . the ratio of total-to-selective extinction, 
Rv_ = A V /E (B — V) = 3.1, and the Galactic extinction curve 
dMisselt et all 199^ ). we can estimate the extinction expected at 10 
eV, i.e., Aio e v and hence rioeiz- For simplicity, we approximate 
the galaxy as a circular disc of radius 5 kpc with uniform surface 
brightness, and the absorber as a poi nt at r = 7.6 kpc along the 
normal to the galactic disc dChen et aPl2005h . to estimate the solid 
angle. Then, for the range of U obtained from the models, the 
expected gas density range is, «h ~30—90 cm -3 . For the above 
density range and /n(h i) obtained from the models, the cloud will 
have size of ~ 0.04—0.5 pc. We have not taken into account the 
angle between the absorber and the normal to the galactic disc 
in this simple calculation, which can cause further attenuation in 


the flux received by the absorber, leading to a lower density and a 
larger cloud size. 

Flence, the absorber can be a halo cloud subject to the radia¬ 
tion field due to a continuously star-forming galaxy, with metallic- 
ity, star formation rate and dust extinction as seen in the candidate 
host-galaxy. Based on the results of our photoionization models 
and our discussion in Section [4] we conclude that the absorber is 
more likely to be tracing gas ejected recently by the galaxy into 
the circumgalactic medium rather than originating in the extended 
disc. Moreover, we note that in all our photoionization model so¬ 
lutions, the iV(Arl) predicted by the model is higher (by ^0.5—1 
dex) than the observed. Argon, being an inert element, is not ex¬ 
pected to deplete on dust, and we find that Ar I is the dominant ion¬ 
ization state of argon in our models. As we pointed out before, Ar I 
depletio n is seen in the FIVCs and in some interstellar sightlines. 
ISofia & Jenkins! 019980 have argued that this could be related to ex¬ 
cess ionization due to hard photons. Our study suggests that for the 
three ionizing radiation fields considered here, Ar depletion cannot 
be explained by a self-consistent ionization model. Ijenkinsl ( l2013h 
suggested that strong Ar depletion could be due to non-equilibrium 
ionization conditions prevailing in the absorbing gas. Alternatively, 
the absorbing gas may be tracing a non-chemically-well-mixed gas 
freshly ejected into the halo from a recently exploded supemovae 
region in the galaxy. In order to make headway on this issue it is 
important to examine several species covering a wide rage of ion¬ 
ization states using future COS observations. 
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Table 4. Summary of the photoionization models for the strong H 2 absorbing component assuming log N(U. 2 ) = 16.61 as discussed in Section [5] 


Models “ 

n H (cm” 3 ) 

/n(h0 

T (K) 

log jV(Cl)/iV(Nal) 

log N(C l*)/N(C 1 ) 

HM12TE 

0.1 

0.5 

160 

2.0 

-2.3 

Galactic Xuv = 0.1—0.5 TE 

[1, 100] 

[0.01, 1.0] 

[20, 50] 

[1.2, 1.5] 

[-1.2,-0.5] 

Galactic XUV = 0.1—0.5 CT 

[3, 30] 

[0.1. 1.0] 

133 

[1.1. 1.2] 

[-1.3,-0.5] 

Starburst TE 

-5.5 6 

0.2 

40 

1.6 

-0.9 

Starburst CT 

[-5.5. -5.0] b 

[0.3, 0.8] 

133 

[1.7, 1.8] 

[-1.1-12] 

Observations 



133^22 

1.8 ± 0.2 

< -0.5 


a TE : thermal equilibrium, CT : constant temperature 
6 log (U) 


5.4 Models with high 1V(H2) 

In Section 13.2.11 we discussed the possibility of a high A r (H 2 ) 
value (10 18,9 cm -2 ) for Component 2. However, in that case the 
molecular fraction will be log f(H 2 ) = —0.57 if we consider all the 
N(H I) to be associated with this component and higher otherwise. 
For the Galactic disc such high molecular fraction generally indi¬ 
cates that the ATHl) is above the threshold (10 20 '' cm -2 ) where 
the H 2 mole cule gets completely self-shielded from interstellar ra¬ 
diation ISavage et alJI 19771) . If we assume the H 2 formation rate 
in this system to be same as in our Galaxy, then the high f(H 2 ) 
would require either a very low photo-dissociation rate (i.e., a much 
weaker radiation field) or a very high density. We remind here that 
the observed T 01 and fine-structure excitation of C does not allow 
the density to be more than few tens of H atoms cm” 3 if the relevant 
parameters are similar to that of the Milky Way (see Section [3.2.2b . 
To consider whether such high molecular fraction is feasible in the 
present system, we ran CLOUDY models with the three background 
radiation fields considered above and stopped the calculations when 
log jV(H 2 ) = 18.9. However, the temperature computed in all three 
cases is much lower (20—30 K) than the inferred T 01 , and hence 
we consider the models with constant temperature. For the extra- 
galactic HM12 background we find that the high N(Hz) can be 
produced for uh = 1 — 100 cm” 3 . However, the IV(Cl) predicted 
by this model is more than 1 dex higher than the observed value. 
For the model with the mean Galactic radiation field, the (V(H 2 ) 
can be produced for hh = 10—30 cm~ 3 and xuv = 0.1. However, 
the 7V(Cl)/(V(Nal) ratio falls ~ 1 dex below the observed value. 
In case of the starburst radiation field, the iV(H I) required by the 
models to produce the ATH 2 ) is higher than the total measured 
N(H I) of the system. Hence, none of the radiation fields considered 
here is able to consistently explain the high N{ H 2 ) solution. This 
supports our argument that the present system is unlikely to host 
such high molecular fraction. However, we emphasize the need for 
higher S/N and higher resolution data to have better handle on the 
models. 


6 RADIO OBSERVATIONS 
6.1 VLBA mas scale imaging 

This sub-DLA is unique in the sense that it is the only known low- 
2 (i.e., z < 1.0) system with a H 2 detection towards a radio-loud 
QSO so that 21-cm absorption observations are possible. Given the 
small inferred sizes from the photoionization models, milliarcsec- 
ond scale VLBI imaging of the background radio source is neces¬ 
sary to measure the covering factor, f c , of the absorbing gas. The 
standard practice, in the absence of VLBI spectroscopy, is to use the 
ratio of the VLBI core flux density to the flux density measured in 


the arcsecond sca l e images to estimate f c (see lKanekar et al .20091 
iGupta et ai1l20 1 2l : [srianand et alJl2012h . We carried out continuum 
observations of the radio source using the VLBA. The spatial res¬ 
olution achieved in our VLBA observations is ~12 mas x 7 mas, 
i.e., ~23 pc x 13 pc at the redshift of the absorber. If the size of 
the H 2 gas is of this order, then the fraction and spatial extent of 
the radio flux density detected in the VLBA image will determine 
the 21-cm absorption detectability. 

The radio_source is identified as a flat-spectrum radio source 
(Healev et al. 2007j), and is unresolve d at arcsecond scales with a 
flux density of 330 mjy at 1.29 GHz jKanekar et al.ll200lh . How¬ 
ever, most of the emission is over-resolved in our VLBA image, 
with only 27% being recovered, assuming the continuum flux den¬ 
sity of the radio source is constant. Our VLBA image (Fig.[5]l shows 
an extended structure with a total flux density of 90 mjy and a 
peak flux density of 59 mjy/beam. Note that we used ‘ROBUST=2’ 
weighting in the AIPS task ‘IMAGR’ to obtain the image shown in 
Fig.0 The largest linear size of the radio source measured from the 
emission detected in the VLBA image is ~70 mas, i.e., ~131 pc at 
the redshift of the absorber. From our discussion on photoioniza¬ 
tion models in the previous section, if the absorbing gas is spher¬ 
ical then the expected extent of the strong H 2 absorber is smaller 
than the extent of the radio emission seen in our VLBA image. If 
we associate the location of the peak emission in the VLBA image 
to the optical source, then the expected covering factor is ^ 0.18. 
However, if we do not assume spherical geometry then we do not 
have any constraint on the extent of the gas in the transverse di¬ 
rection to our line of sight. If we assume the gas to cover all the 
emission seen in the VLBA image, then the covering factor is 0.27. 
In the following section we discussion the implications of this for 
the detectability of 21-cm absorption. 

6.2 21-cm absorption 

In the ATCA radio spectrum j(with a_ channel resolution of 1.8 
kms” 1 ) reported by IlCanekar et al.1 ( bOQlll . no 21-cm absorption 
is seen at the z a b s of Component 2 at the level of 1.56 mjy (the 
spectral rms), while a line is tentatively detected at 3.3<r at z a b s = 
0.10097 in a spectrum smoothed to 9 kms” 1 . The 21-cm optical 
depth t 21 is related to the 7V(Hl) (cm” 2 ) and T S (K) as: ./V(Hl) = 
1.823 x 10 18 (T s // C ) f T 2 idv, where the Hi 21-cm line is as¬ 
sumed to be optically thin. In the galactic ISM. the kinetic temper¬ 
ature^ measured from H 2 usually follows the HI spin temperature 
jRov et al]|2006l ). We assume this to hold for the present system 
as well. If all the 7V(H I) were associated with Component 2, we 
would expect a ~3 mjy line, assuming f c = 0.27 and a line width 
of 5 kms” 1 , which is ruled out at 2cr. Hence, the H 2 gas cannot 
have all the N{ HI) associated with it and cover all the radio emis¬ 
sion in the VLBA image. However, if the H 2 gas were to have ^ 
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Figure 5. Contour plot of the VLBA image of J0441—4313 at 1.4 GHz. 
The rms in the image is 0.15 mJy beam -1 . At the bottom of the image the 
restoring beam is shown as an ellipse. The beam size is 0.012" x 0.007". 
The image centre is at RA = 04 h 41 m 17.3367 s , Dec = -43°13'45.4394". 
The contour levels are plotted as 0.46 X (-2,-1,1,2,4,8,...) mJy beam -1 . 


50% of the total N( HI) or cover only part the VLBA emission, the 
present spectrum would not be sensitive enough to detect 21-cm 
absorption from this gas. A factor 3 higher S/N spectrum will al¬ 
low us to detect or constrain the 21-cm absorption at a significance 
of 3 <t. 

The weak 3.3<r 21-cm detection lies within ~15 kms' 1 of 
the weaker H 2 component, and kms -1 of the weaker compo¬ 
nents of C I, Ca II and Na I. Recall that the measured column den¬ 
sities of metals and dust depletion in this component are not very 
different from that of the strong H 2 component. Usually the ab¬ 
sence/weakness of H 2 and C I ar e asc ribed to high temperatures 
and low densities dSrianand et al.ll2005h . In this case we know that 
the temperature cannot be very high, since the H 2 , C I, Na I absorp¬ 
tion are weak, and the ratio of column densities of the J = 0 & 1 
levels of H 2 in this component suggests that the te mperature is 
100 K. If we consider the optical depth reported bv lKanekar et ali 
d200lh . and assume f c = 0.27 and all of the iV(H I) to be associ¬ 
ated with the weak H 2 component, then T s will be ~90 K. How¬ 
ever, from the observed metal content across the components and 
our photoionization models, the weaker component is not likely to 
have all the 7V(Hl) associated with it. For /n(hi) = 0.5 and f c = 
0.27, the T s will be ~45 K. In case the covering factor is lower than 
0.27, the inferred T s will be even lower. It will be surprising to have 
such a low temperature in a gas with N (H I) ~ few 10 19 cnD 3 . It is 
most likely that one will need higher densities and a lower radiation 
field. In that case it will be interesting to ask why CI lines are weak 
in such a cold and dense gas. The discussions presented in this sec¬ 
tion clearly bring out the importance of having a mu ch high er S/N 
spectr um to confirm the 21-cm detection reported bv lKanekar et all 
fcOOlt) . and have a stronger constraint on the absence of 21-cm ab¬ 
sorption in the strong H 2 component. In addition, better constraints 
on the C I* and C II* absorption will allow us to learn more about 
the physical conditions of the cold gas in this interesting absorber. 


7 SUMMARY 

We have carried out a detailed analysis of the cold gas phase in 
the sub-DLA at 2 a bs = 0.10115 towards 10441—4313, in which H 2 
absorption has been detected bv lMuzahid et alJd2014l) . This unique 
system allows us to study the physical conditions in the absorbing 
gas using both H 2 and 21-cm absorption simultaneously for the first 
time at z < 1. Below we summarize the main results of our study. 

• The H 2 absorption arises from one strong component at z a b s = 
0.10115, with log ATH 2 ) = 16.61 ± 0.05, along with another 
weaker component at z a b s = 0.10091, with log 7V(H2) ^ 15.51 ± 
0.03. We note that the /V(H 2 ) measurement of the strong H 2 com¬ 
ponent is uncertain due the medium resolution of the COS spec¬ 
trum, and that the A r (H 2 ), in principle, can be much higher if the 
actual b parameter is much smaller. However, the absence of H 2 
absorption from higher J levels and of HD absorption, as well 
as our photionization modelling, do not seem to support a high 
iV(H 2 ). The excitation temperature (133^22 K) measured in the 
strong component for the range log ATH 2 ) = 16.6—18.9 is similar 
to what is expected from a CNM phase. 

• The average metallicity of the absorber is found to be twice 
solar and the dust depletion moderate ([Fe/S] ~ —0.49). We do not 
find any significant variation in depletion across the two metal com¬ 
ponents detected in the COS spectrum. The strong H 2 component 
accounts for 55—65% of the metal column densities. Interestingly, 
the stronger components of Cl, Nal and Arl are coincident with 
this component. 

• The sub-DLA is known to be associated with a star-forming 
galaxy at a projected separation of ~ 7.6 kpc, just outside the opti¬ 
cal galactic disk l ichen et ailllOOSl) . We do not find any metallicity 
gradient be tween the host-galaxy and the sub-DLA as suggested by 
IChen et al.1 1 120051) . Moreover, the [N/a] ratio measured in both are 
similar. Hence, the galaxy is likely to either have a metal-rich neu¬ 
tral disk extending to at least ~13 kpc (the galactocentric radius 
of the sub-DLA) or have undergone recent periods of metal-rich 
outflow. 

• From photoionization modelling, we show that the observed 
column densities of H 2 , Cl and Nal in the z a b s = 0.10115 com¬ 
ponent can be consistently explained using a radiation field due to 
a continuously star-forming galaxy, with metallicity, star formation 
rate and dust extinction as measured in the associated galaxy. Al¬ 
ternatively, if the absorber is tracing the extended galactic disc then 
the radiation field has to be weaker than half the mean Galactic radi¬ 
ation field. However, the models suggest that the absorber is more 
likely to be tracing gas in the galactic halo than in the extended 
galactic disc. We note that, the measured column densities are also 
consistent with the extragalactic background radiation if the gas 
has 71 h ~ 0.1 cm -3 and L ~70 pc. However, from the presence of 
a star-forming galaxy at ~7.6 kpc, we argue that this absorber is 
more likely to be a halo cloud. In that case, using simple approxi¬ 
mations we obtain for this gas cloud, uh ~ 30—90 cm~ 3 and L ~ 
0.05—0.4 pc. 

• The measured Ar I absorption strength in this system is much 
weaker than what is expected if there is no depletion or ionization 
effects, and than what is predicted by all our photoionization mod¬ 
els. Similar depletion is seen in case of HVCs in our Galaxy. The 
possible explanations could be non-equilibrium ionization condi¬ 
tions in the absorber or the gas not being chemically well-mixed 
due to recent ejection from supernovae region in the associated 
galaxy. 

• This sub-DLA presents a fortuitous case of a H 2 absorber to¬ 
wards a radio-loud QSO which can be searched for 21-cm absorp- 
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tion. We present the VLBA mas scale image of the background 
radio source, which shows it to be over-resolved, with only 27% of 
the arcsecond flux being recovered. The radio emission has a total 
flux density of 90 mJy, a peak flux density of 59 mJy/beam and an 
extent o f ~131 p c at the redshift of the absorber. 

• iKanekar et alj < 2001 1 report a tentative detection of weak 
21-cm absorption that arises within ~15 kms -1 of the weaker 
H 2 component. However, 21-cm absorption from the stronger H 2 
component is absent at the level of 2cr, indicating that either the 
7V(H I) associated with this component is ^ 50% of the total mea¬ 
sured iV(H I) or that the covering factor of this gas is < 0.27, con¬ 
sistent with results from our photoionization models. 

• The reported 21-cm absorption from the weaker H 2 compo¬ 
nent indicates that T s ^ 90 K in this gas. The actual value could be 
much lower if a good fraction of the observed N( H I) is associated 
with the strong H 2 component. However, the weakness of H 2 , Cl 
and Na I absorption from such cold gas is puzzling. A higher S/N 
21-cm absorption spectrum and a UV spectrum with better cov¬ 
erage of species like Cl, Cl* and C II* are essential to put more 
stringent constraints on the conditions in this absorbing gas. 

Such detailed analysis of similar individual systems is essential 
to further our understanding of the properties of cold gas present 
around galaxies, and in particular, of the connection between 
H 2 and 21-cm absorption. The results presented here will also 
facilitate interpretation of the results of large surveys of H 2 and 
21-cm absorption. 
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